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Abstract. Recent electronic structure calculations for the title compounds 
performed by Wu et. al. [nj are critically reconsidered, applying high precision full- 
potential bandstructure methods. It is shown that the bandstructure calculations 
presented by the authors contain several important inconsistencies, which make their 
main conclusions highly questionable. 



In a recent paper Wu et. al. [Q presented bandstructure calculations for the quasi one- 
dimensional CuO-chain compound SrCuC>2 and the quasi two-dimensional material 
CaCuC>2, both being of prototypical character and therefore of general interest. Wu 
et. al. used a full-potential linear combination of atomic orbitals method g] in the 
framework of the local spin density approximation (LSDA) and included on-site 
Coulomb interaction corrections (LSDA+t/). The authors of Ref. [[t] claim that on the 
basis of their full-potential band structure experimental findings can be well fit with 
an U of 5 eV, significantly smaller than U values reported in previous calculations || . 

However, there are obvious inconsistencies and important differences between 
the calculations of Ref. |l| and previous studies ||, [|, [j| ||, concerning (i) the 
proper symmetry in fc-space, (ii) the widths and the orbital character of the shown 
bands, (iii) the total (DOS) as well as the partial densities of states (PDOS). 
Therefore, we reinvestigated the electronic structures of CaCu02 and SrCuC-2 using 
two independent, well basis converged full-potential bandstructure methods to find 
out whether or not the differences mentioned above could be understood as a 
consequence of the differences between a full-potential [jlj and the earlier non-full- 
potential calculations ||, |[ [|. We carried out LSDA bandstructure calculations for 
CaCu02 within a full-potential minimum-basis local-orbital scheme (FPLO) |7| and 
within a full-potential linearized augmented plane wave (FLAPW) scheme both 
in scalar relativistic versions. (We note that relativistic effects are in the order of 0.1 
eV only.) In the FPLO-scheme, modified Ca 3d, 4s, 4p, (Sr 5s, bp ,4d), Cu 3d, 4s, Ap, 
and O 2s, 2p, 3d states were used as valence states for CaCu02 (SrCuC^), the lower 
lying states were treated as core states. The WIEN97-code || employs local orbitals 
(LO) to relax linearisation errors and to treat the 0-2s and semicore Cu-3p and Ca-3s, 
3p states. Well converged basis sets of over 500 APW functions plus LOs were used. 
The radii of the atomic spheres in the latter case were 1.8 a.u. for all atoms. The 
basic calculations were performed with 125 and 90 fc-points in the FPLO-scheme and 
in the WIEN97-code, respectively, for the irreducible part of the Brillouin zone using 
the tetrahedron method. We itasize that the numerical convergence (with respect 
to the number of fc-points Nk, the valence basis set, the potential and the density 
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Figure 1. Total DOS (a,b) and partial DOS (c,d) of CaCu02 calculated within 
various bandstructure schemes (see text). 



representataion) of all calculated properties was carefully checked. Following Ref. 
we will discuss first CaCu02 and afterwards SrCuC>2. 

CaCuC>2 First, we will concentrate on the bandstructure and then on the DOS. 
Both bandstructures obtained with the FPLO and WIEN97 codes agree excellently 
with each other (see Fig. 1 (a)) and with previously published results [j| . Considering 
the bands in Ref. |l[] one realizes the following points: The authors show a 
bandstructure with orthorhombic symmetry for the tetragonal crystal structure || 
(note the different dispersions in Fig. 1(a) of Ref. Q along the r-(100) and r-(010) 
direction, respectively, which must be equivalent for the tetragonal case). Also, the 
band degeneracies at symmetry points are incorrect. As a consequence, the number 
of degeneracies in r-(OOl) direction (c-direction) differs from all other calculations 
(11 different bands instead of 8 different bands allowed by the crystal symmetry). 
Contrary to our results and to the results of Refs. (J, |) , the authors find an additional, 
third band with sizable dispersion with its maximum at (001). We analyzed the 
orbital character of our bands, in particular to find out which states are responsible 
for the relatively large dispersion in c-direction of about 1-2 eV discussed also by 
Mattheiss et al. ||. Our calculations show that these two strongly dispersive bands 
have predominant O 2p z character with a small admixture of Cu 3d 3z 2_ r 2 states. In 
contrast, the c-dispersion of the antibonding band essentially made up by O 2p-a 
and O 2p-n orbitals with the Cu 3d x 2_ y 2 orbital is only about 350 meV. Just these 
states mediate the magnetic coupling between different layers. This confirms the quasi 
two-dimensional character of the magnetic Hamiltonian. 

In Fig. 1(b) we compare our results for the DOS with those of Fig. 2(a) in Ref. Q. 
The width of their pd-complex is too large by about 2 eV. The authors attribute this 
discrepancy to their choice of an ionic basis, which only means that their calculations 
are not basis set converged. The reason for the large discrepancy between the DOS of 
Ref. |l) and our's is evident from Figs. 1(c) and 1(d). Due to the ionic orbital basis used 
in Ref. Q , the O 2p states are shifted downwards by about 2 eV and the hybridisation 
with Cu 3d states is consequently reduced. We attribute simply the reduction of the 
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Figure 2. FPLO total DOS for SrCuC>2 of the pd-complex (left panel) and the 
zoomed region near the Fermi level for different number of fc-points (right panel), 
calculated in the doubled unit cell. 



bandwidth in the LSDA+C/ in the calculation of Ref. 0] to a downwards shift of their 
Cu 3d states towards the incorrectly positioned oxygen 2p states. 

SrCuC>2 For our paramagnetic calculation the resulting DOS is in excellent 
agreement with the DOS reported in Refs. |6, 11 1. One should note that the occurrence 
of van-Hove singularities at the band edges of the antibonding band, due to the 
nearly one-dimensional electronic structure of the compound, depends critically on 
the sufficiently large Nk used in the calculation. 

For the supercell calculation, we doubled the unit cell along the chain direction 
and started the self-consistent calculation with an antiferromagnetic arrangement of 
the Cu spins along the chains Q . To describe properly some peculiarities related to 
the nearly one-dimensional electronic structure, we made several calculations varying 
Nk- In particular, we enlarged Nk along the chain direction. The results are shown in 
Fig. 2. Due to the nearly ideal one-dimensional dispersion of SrCu02, the calculation 
results in an erroneous insulating groundstate for an insufficient Nk- At least for 
Nk > 250, the artificial gap and the related singularities disappear and the results 
converge towards those of our paramagnetic calculation. Possibly, the gap of 0.55 eV 
in Ref. jlj can be attributed qualitatively to an insufficient number of fc-points. The 
reported relatively big magnetic moment of 0.33 [Ib is related to this artificial gap 
and to the unusualy small hybridisation of the Cu 3d states with the O 2p states. An 
antiferromagnetic solution has been reported also by other authors Q though with an 
even much smaller gap and extremely small magnetic moments for the SrCu02-system. 
Again, we attribute this gap (notably smaller than that in Ref. Q due to a larger Nk) 
to a still too small number of fc-points. To make this point more clear let us consider 
schematically the band structure in the folded zone picture in T - (001) direction near 
the antiferromagnetic Bragg point (see Fig. 3). The weak hybridisation between the 
two subchains of the double chain leads to two slightly split bands. After the folding 
of these two bands in the doubled unit cell, new crossing points appear close to the 
new symmetry plane. Therefore, DOS routines result in a wrong interpolation if the 
crossing points are not at calculated fc-points. The mentioned above splitting A max 
« 150 meV (see Fig. 3) provides an upper bound for the artificial gap |L3|. Hence, for 
supercell calculations, especially for quasi one-dimensional electronic structures, one 
has to be very careful choosing Nk in applying standard interpolation methods. 

In Fig. 4(a) of Ref. g) one finds the Cu 3d states in the PDOS with a too small 
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Figure 3. Schematic picture of the folded antibonding bands in T - (001) direction in 
SrCu02- The solid lines are the correct bands, the dashed lines show the interpolated 
bands for the k points given by the open circles. For explanation see text. 

width as already discussed for CaCu02- Moreover, for an orbital projected DOS, it is 
reasonable to use symmetry related orbitals with their quantization axis perpendicular 
to the Cu04 plaquette. In this representation only the Cu 3d x 2_ y 2 orbital contributes 
considerably to the antibonding band (instead of two Cu 3d orbitals shown in Fig. 4 
of Rcf. @), as found for CaCu0 2 (Fig. 1(a) of Ref. 0). The LSDA+J7 procedure 
used in Ref. [[[] depends on the basis set representation and would require therefore 
the application of the same local orbital symmetry for both compounds in order to 
ensure a proper comparison between them. 

To summarize, we discussed the main differences between the results of Wu et 
al. and our/or previously published data || ^, ||, [| basic paramagnetic LDA 
calculations contain several serious inconsistencies, all conclusions reported in Ref. |jj 
with respect to LSDA+C/ are highly questionable. In particular, the large gap in 
SrCu02 at moderate U might result from the artificial gap found in their LSDA 
calculations. 
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